Inverse velocity statistics in two dimensional turbulence 
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We present a numerical study of two-dimensional turbulent flows in the enstrophy cascade regime, 
with different large-scale forcings and energy sinks. In particular, we study the statistics of more- 
than-differentiable velocity fluctuations by means of two recently introduced sets of statistical es- 
timators, namely inverse statistics and second order differences. We show that the 2D turbulent 
velocity field, u, cannot be simply characterized by its spectrum behavior, E(k) cx k~ a . There 
exists a whole set of exponents associated to the non-trivial smooth fluctuations of the velocity field 
at all scales. We also present a numerical investigation of the temporal properties of u measured in 
different spatial locations. 
PACS: 47.27.Eq, 47.27.Gs 
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I. INTRODUCTION 



Many natural phenomena display complex fluctuations 
over a wide range of spatial and temporal scales. Com- 
plexity usually manifests in the non-Gaussian properties 
of probability distribution functions (PDF). When PDFs 
at different scales do not collapse by a simple rescaling 
procedure one speaks about intermittency jl ] . Such non- 
trivial rescaling properties may be exhibited by PDFs' 
tails or peaks, or both When intermittency manifests 
in the PDF's tails, it means that regions of very intense 
bursting activity are present. This is typical of three 
dimensional turbulent flows, where the velocity field is 
strongly intermittent and rough jy . 
However, there are examples of other important natural 
phenomena which develop simple PDF's tails but non- 
trivial PDF's cores. PDF's peaks are associated to lam- 
inar fluctuations, i.e., "smooth" variations of the field. 
A physically relevant example is offered by two dimen- 
sional turbulent flows where the presence of long living 
coherent structures, e.g., vortices, is very well known (see 
Figure 1). Two dimensional turbulence is characterized 
by two different transport processes : an inverse energy 
cascade from the forcing scale to larger scales and a di- 
rect enstrophy cascade from the forcing scale to smaller 
scales ||. Inverse energy cascade shows a non intermit- 
tent Kolmogorov 1941 scaling for the velocity field j| [|. 
On the contrary, in the direct cascade non-trivial vor- 
ticity fluctuations have been observed in dependence on 
the large scales characteristics of the flow [0-|o| . In addi- 
tion, velocity fluctuations in the direct enstrophy cascade 
regime are particularly interesting for geophysical and as- 
trophysical sciences In this regime, the velocity field 
is diffcrentiable, therefore the standard analysis (custom- 
arily applied in 3d turbulence), based on moments of ve- 
locity increments (the so-called structure functions) is 



poorly informative. Indeed structure functions are dom- 
inated by the differential component of the signal: 



S p (r) = ([s(x + r) - s{x)] p ) ~ r F 



(1) 



or the u v ve- 



where with s we indicate either the u x 
locity fields component. It is worth stressing that the 
scaling behavior (|l|) does not imply that the velocity 
statistics is trivial. For example, it is well known that 
in the enstrophy cascade regime the energy spectrum 
shows a power law E(k) cx k~ a with a > 3, which is 
the signature of significant more-than-differentiable ve- 
locity fluctuations. Hence, subdominat contributions to 
the s{x + r) — s(x) cx r behavior must be present and, in 
principle, detectable. The triviality of the scaling (|l|) it 
is just the consequence of not having chosen the suitable 
observable. Therefore, to extract interesting information 
on the statistics of smooth signals, new statistical tools 
are needed. 

Recent contributions have shown that laminar events are 
optimally characterized in terms of their exit-distance 
statistics, also known as inverse statistics |i^-|l6[]. In 
a nutshell, in such approach one "inverts" the usual 
way of looking at signals. Standard analysis studies the 
statistics of signal increments over a certain spatial (or 
temporal) interval; the exit-distance approach looks at 
the statistics of spatial (temporal) intervals necessary 
to observe a given signal increment. Another possibil- 
ity to study smooth signals is to eliminate the differen- 
tiable contribution by looking at signal Second Differ- 
ences (SD), i.e., (s(x + r) — 2s(x) +s(x — r)) as suggested 
in P. 

In this paper we extend a previous exploratory investiga- 
tion |6| of the inverse statistics of velocity fields in the 
enstrophy cascade regime of 2D turbulence, and we com- 
pare it with results obtained by Second Difference statis- 
tics on the same flows. We present both exact analytical 
results for the exit-distance probability density functions 
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of ID Gaussian signal, and a set of numerical investi- 
gations of spatial and temporal statistics of 2D turbu- 
lent flows. The main result is the identification of highly 
non-trivial contributions to the more-than-differentiable 
velocity fluctuations. We also introduce a set of expo- 
nents which characterizes smooth behaviors beyond that 
of the energy spectrum, a, E(k) oc k~ a . 
The paper is organized as follows. In section [ij], we recall 
some known results on 2D turbulent flows in the p rese nce 
of a drag mechanism at large scales. In section III, we 



introduce the main observables, i.e., the inverse structure 
functions and the the Second Difference structure func- 
tions : we first apply them to the analysis of stochastic 
signals with a given spectrum E(k) ~ k~ a , for which we 
are able to establish some exact results. Then in sec- 
tion IV , we present the spatial statistics of laminar fluc- 



tuations of the two dimensional velocity field u obtained 
by direct numerical simulations (DNS). In section [y], we 
perform a temporal analysis of the velocity field on fixed 
spatial locations. Section VI is devoted to conclusive re- 
marks. 




FIG. 1. A snapshot of the vorticity field. Colors are coded 
according to the intensity of the vorticity field from the min- 
ima of uj (black) to the maxima (clear yellow,white). DNS 
have been performed by a standard dealiased pseudo-spectral 
algorithm, over a double periodic square domain of size 
L — 2ty, at resolution 512 2 and 1024 2 . As customary en- 
strophy is dissipated at small scales with an hyper-viscosity 
of order 4, while energy is removed at large scales, to avoid pil- 
ing up on the smallest mode, using different drags (see table). 
We considered a Gaussian, white-in-time large-scale forcing 
restricted on wave numbers, 4 < \kf \ < 6. 



II. TWO DIMENSIONAL TURBULENCE 

As far as the inertial range of scales for the enstro- 
phy cascade of two-dimensional turbulence is concerned, 
previous experimental, theoretical and numerical studies 



have shown that the statistics is strongly influenced by 
large-scale phenomena. Indeed more than smooth spectra 
E{k) ~ k~ a with a > 3, depending on the characteris- 
tics of the forcing and of the large-scale dissipation, have 
been reported Recently, new results have clarified 

the problem in the special case of the large scale energy 
sink given by a linear (Eckman) friction p|]lO|]. We re- 
call that the presence of an energy sink at large scales 
is conceptually justified by the necessity of avoiding the 
pile up of energy on the gravest mode as a result of the 
inverse energy cascade Q and it is physically motivated 
in terms of the friction to which a fluid is subjected in 
the Eckman layer ■ 

The strong influence of large-scale phenomena in the 
whole enstrophy cascade range is believed to be a conse- 
quence of non local interactions (in Fourier space). An- 
other property associated to the enstrophy cascade is the 
velocity field smoothness. The aim of this paper is to 
discuss a new set of observable suitable to highlight the 
statistics of all those fluctuations which appear as a sub- 
leading contribution to the smooth differentiable behav- 
ior, u(x + r) — u(x) ~ r. 

Let us now briefly fix the notation. In terms of the scalar 
vorticity u> — V x u, the equation of motion can be writ- 
ten as 



dtOJ + u ■ Vw 



v q A q uj 



P P A-Pu} + F, 



(2) 



where v q and f3 p (q, p > 0) are the coefficients of the gen- 
eralized dissipations, namely the hyper-dissipative and 
the hypo-friction terms respectively. The former removes 
enstrophy at small scales and the latter removes energy 
at large scales. F is the vorticity source term acting at 
large scales. In Fig. 1, we show a typical snapshot of the 
vorticity field obtained by direct numerical simulation of 
Eq. (j^). As one can see the vorticity field is characterized 
by filamental structure over a wide range of scales. 
According to the classical prediction || , the velocity field 
should exhibit a Batchelor-Kraichnan spectrum, E(k) ~ 
fc _3 (ln(fc))" 1 / 3 . The dimensional estimate is observed 
in a bunch of numerical and experimental measurements 
However, in the literature there are reported nu- 
merous situations [@-|l0| where different velocity spectra 
have been measured : E{k) ~ k~ a , with the exponent 
a larger than 3 and dependent on the forcing and drag 
mechanisms. In the case of linear friction, (p = 0), it is 
known that vorticity statistics is intermittent. In such a 
case, it has been recently clarified ||[l0| that, at scales 
small enough, vorticity behaves as a passive scalar. In 
addition, the dependency of the spectrum slope on the 
linear-friction coefficient has been understood . Ex- 
cept for the situation with a large-scale linear friction, 
there is no general theory for the scaling properties of 
2D turbulent flows in the presence of different large-scale 
drag mechanisms (see also p2[). 

Let us therefore present the way we analyzed 2D turbu- 
lent flows with general large-scale physics and the related 
results on their statistics. 
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III. INVERSE AND DIRECT STATISTICS FOR 
SMOOTH SIGNALS 

In this section we introduce the inverse statistics and 
the second difference structure functions. We start ap- 
plying them to the analysis of stochastic one-dimensional 
signals with a given spectrum E{k) ~ k~ a . For the 
sake of simplicity we limit our discussion to signals with 
3 < a < 5, for which we are able to establish some exact 
results. To be more precise, we consider smooth random 
signals built as follows 



T ( p\Ss) = (r p {Ss)) ~ 8s x{p) 



(4) 



i(kx+e k ) 



(3) 



where \s(k)\ 2 ~ k~ a and 9k are random phases, uni- 
formly distributed in [0, 2ir). When 3 < a < 5 the signal 
is smooth but only one time differentiable. Hence, mo- 
ments of its differences over any increment r always pos- 
sess a differentiable scaling (|l|), while moments of order 
p < — 1 do not exist. 



A. Inverse statistics 

For a generic smooth one-dimensional signal s(x), look- 
ing at inverse statistics consists in measuring moments of 
the distance, r(5s), necessary to observe in the signal a 
double exit (forward and backward) through a barrier 8s. 



S(X)A 




FIG. 2. A pictorial representation of the exit-distance 



where the average is taken with respect to the random 
choice of the point xq |23|] . For the case of simple signals 
such as (||), a rigorous estimate of the scaling exponents 
of inverse statistics moments can be derived as follows. 
If the signal spectrum is E(k) ~ k~ a with 3 < a < 5, 
we can write the signal increment as 



s(x + r) — s(x) 



ds(x) 
dx 



- c(x) r h 



(5) 



Here we have only kept the two most important scal- 
ing behaviors: 0(r) because of the differentiability and 
0(r h ) from the spectrum exponent. The scaling expo- 
nent 1 < h < 2 is related to the spectrum slope by the 
dimensional relation a = 1 + 2h, while c{x) is a continu- 
ous function of x. By studying the exit event, in the limit 
of a small barrier height, we may observe two different 
kinds of event. The first, with probability one, is the 
differentiable scaling r(Ss) ~ 5s. The second, observed 
at those points where the first derivative vanishes, is the 
subleading behavior, 0(r h ), in (||). 

One may estimate the probability of this second situation 
as follows. With 3 < a < 5 the first derivative is a self- 
affine signal with Holder exponent £ = (h— 1) < 1, which 
vanishes on a fractal set of dimension D = 1— £ = 2 — h. 
Therefore, the probability to see the sub-leading term 
0(r h ) dominating the exit event in (^|) is given by the 
probability to pick at random a point on a fractal set of 
dimension D, i.e., 



P(r ~ {Ss) 1 ^) - r x ~ D = {Ss) 1 ' 1 / 11 . 



(6) 



Taking into account both events, we end with the follow- 
ing bifractal prediction for inverse statistics moments : 



P 



1 - - 



(7) 



From (Q) , we conclude that laminar differentiable fluctu- 
ations influence the inverse statistics only up to moments 
of order p = 1; for larger p, the PDF is dominated by the 



method. X; and Xj are two points picked at random and sub-dominant behavior, (s(x + r) — s(x)) 



ri and rj are the corresponding exit distance from the barrier 
8s. 

We fix a value for the signal fluctuation, 8s, then we 
pick at random a reference point xq and measure the 
first forward (|s(xo + rr) — s(x a )\ > 8s) and backward 
(|s(xo — Tb) — s(xo)\ > 8s) exit from the barrier, r(8s). 
Then we put r(x , Ss) = + r/. See Fig. 2 for a picto- 
rial view of the method. Repeating the observations for 
many point xo and for different barrier heights, we can 
define the inverse structure functions Q as 



\ In other 

words, the extrema of the signal play the role of sin- 
gularities for the inverse statistics : close to the extrema, 
events with much longer exit distances are observed when 
Ss —>■ 0. For one-dimensional signals as (||) the prediction 
is verified with high accuracy (see Fig. 3). 
In the general multiaffine case, signal increments scale as 
S r s(x) ~ r h ^ with probability P r (h) ~ r i-£ ) ('i) j where 
the function D(h) can be interpreted as the fractal di- 
mension of the set where the Holder exponent h is ob- 
served Q. 

For such a signal, it is possible to obtain fl4[13| a link 
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between the inverse statistics exponents, x{p)^ an d the 
fractal dimension, D{h): 

. f p+l-D(h) \ 

^p) = i T{ h )■ (8) 

In the case of the smooth signal ((H) , one can see that (||) 
coincides with the bifractal prediction ([?]) , as soon as we 
write D(h) = 2-h for h = (1, (a - l)/2). 
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FIG. 3. Scaling exponents x(p) f° r the ID signal ([]) with 
a — 4. The solid line gives the bifractal behavior. Moments 
have been computed using 10 3 realizations of the signal ([]) 
with 2 17 modes; for each realization 2 12 starting points, xo, 
have been taken at random. 



B. Second Difference structure functions 

Another way to eliminate the trivial differential scal- 
ing and extract some statistical information from smooth 
signals has been suggested in jl7). The basic idea is 
to consider moments of the second difference A,.s = 
(s(x + t) + s(x — r) — 2 s(x)), so that the differentiable 
contribution, 6 r s cx r, is automatically eliminated. For 
the signals under investigation, we have that at the lead- 
ing order A r s ~ r h with 1 < h < 2, and moments behave 
as 

4^(r) = <|A rS | p >~r*». (9) 

In the monofractal case (globally self-similar signals), one 
expects z p = ph. The analysis done for the same stochas- 
tic ID signal of (||) with h = 1.5, confirms this expecta- 
tion (see Fig. 4). 

In the general case, i.e., when many more-than- 
differentiable fluctuations are present, the scaling expo- 
nents z p are non-trivially related to the distribution of 
the h exponents. The difficulty to give a multifractal 
prediction for A r s increments stems from the fact that it 
is a three-point quantity, depending on the simultaneous 
fluctuations between (x, x — r) and (x, x + r). There- 
fore, to draw the multifractal picture, we would need in 
addition a complete control of the spatial correlations. 



Something which is in principle feasible |2^,|26| but which 
goes beyond the aim of this paper. 
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1 2 3 4 5 6 



FIG. 4. Scaling exponent for the ID signal (g) with the 
same parameters as in Fig. 3, the straight line shows the ex- 
pected behavior z v = hp with h = 1.5. The inset shows the 
local slope for p = 4. 



IV. SPATIAL STATISTICS IN 2D SMOOTH 
VELOCITY FIELDS 



Let us now analyse the inverse and second difference 
statistics of the two dimensional velocity field obtained 
by DNS of the Navier-Stokes equation (||). We performed 
four different sets of numerical experiments, with peri- 
odic boundary conditions on a spatial grid of 1024 2 col- 
location points. In all of them, we considered a Gaussian 
forcing, 5-correlated in time, with support in a restricted 
band of wave numbers 4 < kf < 6. 



LU 




FIG. 5. Log-log plot of the velocity spectra for two differ- 
ent drag coefficients with the Eckman linear friction, run B 
(middle) and run C (bottom), and with a hypodiffusive fric- 
tion, run D (top). Straight-lines correspond to the best fit 
power laws, A; -3,38 , k~ 3 ' 74 and fc -3 ' 26 respectively. Run A is 
not shown because it is almost indistinguishable from run D. 
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In three out of four simulations, we used Eckman lin- 
ear friction, i.e. p = in (|^) with different coefficients 
(simulations A,B,C, in the following). In the fourth run, 
we used a hypo-diffusive term at large scales, p = 2, (re- 
ferred as case D in what follows). Table 1 is a summary 
of the DNS parameters, together with the best-fit spec- 
trum exponent a for all runs. 

In Fig. 5 we show the averaged velocity spectrum for run 
B,C and D (run A gives a slope almost coincident with 
that of run D). By comparing them, it is evident that 
the spectrum slope depends on both the drag coefficient 
(runs A,B and C) and on the drag mechanism (run D). 
Evidently, we are in presence of large-scales effects which 
somehow affects small scales velocity fluctuations. Let 
us try to quantify this statement by using the inverse 
statistics analysis. First we compare the inverse struc- 
ture functions measured on several snapshots of the DNS, 
with those obtained after randomization of all velocity 
phases on the same frames. The rationale for this test 
is to investigate the importance of correlations between 
fluctuations at different wave-numbers and therefore the 
"information" content brought by coherent structures in 
2D turbulent flows. 

If we look at a one-dimensional cut of the velocity field, 
before and after phases randomization, it is rather diffi- 
cult to distinguish the original DNS field from that-one 
with randomized phases. This is due to the steepness 
of the spectrum, i.e., only few modes dominate the real- 
space configuration. Despite the apparent similarity big 
differences arise when looking at inverse moments. 
Because of the limited numerical resolution, the only 
quantitative statements one can give are for relative scal- 
ing properties. Therefore, we measure scaling laws of the 
inverse statistics by plotting all moments T^ p \5u) versus 
a reference one, s ay T^idu). This is the same tech- 
nique called ESS p7fl , fruitfully applied in the analysis 
of 3D turbulent data with the aim of re-absorbing some 
finite size effects and extracting scaling information also 
at moderate resolution. Therefore, we concentrate on the 
following relative scaling properties : 

T (p) (Su) oc (T^(<kt)) x(p)/x(2) ■ 

In Figs. 6 we summarize our findings. Inverse moment 
exponents, x(p)/x(2), measured on the turbulent fields 
with randomized phases follow the bifractal prediction 
(0) , with the value of h extracted from the spectrum (see 
Table I). Conversely, the longitudinal and transversal 
inverse-statistics moments without phases randomization 
show anomalous scaling laws, which deviate from the 
bifractal law given in (In). In Figs. 6, we show the curve 
x(p)/x(2) for both randomized and non-randomized 
transversal exit moments for run C and D. For p < 1, 
the statistics of the randomized data and that of the 
turbulent data almost coincide being those moments 
(with < p < 1) dominated by the laminar fluctuations 
u(x + r) — u(x) ~ r. To better appreciate differences in 
the scaling curves, we show in the inset of Figs. 6 the 



' (Su) , for the random- 



local slopes of T^ 4 ) (Su) versus 

T (2) 

ized and non randomized data. 

The following scenario can be drawn. Randomized data 
follow the bifractal prediction, while the non randomized 
ones are definitely different and display anomalous scal- 
ing. Moreover, the anomalous scaling is present for all 
choices of the drag mechanism. 
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FIG. 6. x(p)/x(2) for exit moments obtained in run C 
(top) and D (bottom) (o), compared with the same data af- 
ter phases randomization (x). In both figures, the solid line 
corresponds to the bifractal prediction. The values of h used 
for the theoretical prediction are obtained from the spectrum 
scaling exponents (see Table 1), namely h = 1.38 (run C) 
and h = 1 .13 (run D). Errors on scaling exponents have been 
computed according to the fluctuations of the local slopes. 
In each inset, it is shown the local slope of T^'(Su) versus 
T^ 2 '(Su) for the true signal (o) and the randomized one (x). 

For Second Difference statistics analogous results have 
been found, that is a monofractal behavior for the ran- 
domized field and an anomalous behavior for the turbu- 
lent one. In Figs. 7 we show the scaling exponents z p for 
run (C) and run (D). Longitudinal and transversal com- 
ponents, within the errors, coincide. The SD analysis 
confirms that the statistics of laminar events for the 2D 
turbulent velocity field displays a complex, multifractal 
structure. 

Concerning the case of run C, i.e., with linear friction, it 
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is interesting to compare the results of the Second Dif- 
ference moments with some recent analytical results [ p8| . 
In ps| ], it is argued that in presence of linear Eckman 
friction, the second and third order (standard) structure 



functions behave as 5*2 (r) 
and S 3 (r) = ((6u(r)) 3 ) ~ 



{Su 2 (r)) 



dr A 



- ar- + br 2+( - a -^ 
being a > 3 the 



spectrum slope, and a, b, d, e some constants. From these 
results, it is easy to extract the exponents of the SD 
moments, i.e., Zi = a — 1 and Z3 = a. Actually our data 
give slightly larger values. Such discrepancies may be 
due to strong finite Reynolds effects which, in 2D, are 
particularly severe due to the interplay between inverse 
cascade and friction in the low k region of the spectrum. 





FIG. 7. (top) Second Difference exponents, z p , for mo- 
ments obtained in run (C) (o) and after randomization (x); 
(bottom) the same but for run D. The straight lines corre- 
spond to the monofractal behaviour z v — hp with h — 1.55 
(run C) and h = 1.20 (run D). 



V. TEMPORAL STATISTICS 

As it is well known, in 3D turbulence we can recast 
the temporal behavior of the velocity field into the spa- 
tial domain via the Taylor hypothesis (frozen turbulence 
hypothesis): the effect of large scales is just that of a 
uniform sweeping which does not modify the small scale 
structures and their energy content. In 2D the absence of 



a time hierarchy rules out such a possibility jj 21|. This 
is also evident by looking at snapshots from numerical 
simulations, which show that the time evolution of the 
dynamics is dominated by stable, long-lived structures 
(see Fig. 1). For such a reason it is non trivial to predict 
the velocity temporal statistics collected in a fixed spa- 
tial location. 

We performed a DNS of @ taking as a large-scale forcing 
a function F of constant amplitude at some characteris- 
tic wave- numbers 4 < \kr\ < 6 and time-independent. 
We performed a long time integration of the 2D N- 
S equations, at resolution 512 2 and collected statistics 
for hundreds large eddy turn over times, estimated as 
teddy ~ 1 / Wmj (details on the numerical simulation can 
be found in f29(] ). Once the system reached a station- 
ary state, we started to collect the time evolution of the 
velocity fields at some specific spatial locations with a 
sampling time T samp ~ 2.5 • I0~ 2 t eddy . 
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FIG. 8. (top) Time recording of the x-component of the ve- 
locity field of the Pi n probe (dotted line) and the p out probe 
(solid line); (bottom) log- log plot of the frequency spectra of 
the following signals: (a) E Piri (zu) of the probe near a coher- 
ent structure; (b) E Pout (vj) of the laminar one; (c) Et f (ro) of 
the time evolution at a particular Fourier mode k, in forced 
wave- number band |fe/| = (4; 6]. 

Some observations are noteworthy. The first one con- 
cerns the ergodicity of the velocity field u(x,t). Tempo- 
ral signals collected at different spatial locations possess 
different probability distribution functions. In particular 
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the range of variations of the local rms values u rms (xQ, t) 
is so wide that we can not average time histories recorded 
at different points. It is difficult to say if waiting long 
enough one would recover, as expected, some stable er- 
godic properties. Certainly, with our statistics we feel 
confident to report results only on local averages, avoid- 
ing to mix temporal evolutions in different spatial loca- 
tions. 

In particular, we chose to report those describing two 
typical spatial situations: one, Pi n , situated in the core 
of a vortical structure, the other, p out , in a laminar re- 
gion. This means that notwithstanding the turbulent 
evolution of the field, the motion of the vortices is so 
slow that probes almost maintain their respective "char- 
acter" (in and out of a vortex) all the simulation long. 
Looking at a sample of the time series recorded by the 
two probes, the signals seem very different and change 
when the probe passes from a laminar region to a vorti- 
cal one (see Fig. 8 (top)). To have a better understand- 
ing, it is useful to consider the frequency spectra E(tu) of 
the signals, calculated from the temporal Fourier trans- 
form of the stationary time correlation function, e.g., of 
C(u x (p),r) = (u x (p,t + T)u x (p,t)}. 
In Fig. 8 (bottom), we can observe the spectra of the 
two probes. At variance from the spatial spectra, it is 
not possible to extract a clear scaling behavior. One can 
only identifies an exponential decay, and a peak region 
located at the frequency zu^ « 0.51/t e ddy (here and in 
the sequel "(I*)" stands for "large scale".) It is easy to 
recognize that T^ L \ defined as = l/m {L \ is the typ- 
ical time scale associated to the large-scale structures, ei- 
ther estimating it from the v orticity content of the largest 
structures T^ L ' = 1/ yj (lu 2 ) or from their typical revolu- 
tion time. In other words, Fig. 8 (bottom) tell us that 
in each spatial point the time evolution is governed by 
the typical oscillation frequency of the forced large-scale 
structures. 

This is confirmed by the comparison of the spectra probes 
with the spectra built from the time correlation of the 
Fourier transformed velocity field u(k,t), at a given 
mode, kf, belonging to the forced wavenumber band. 
Indeed, all spectra posses a peak at frequency w^ L \ 
We now pass to investigate direct and inverse statistics 
of u(t). Direct structure functions behave trivially for 
both probes, S p (t) = ([u(x , t + r) — u(x , t)] p ~ c p t p , 
where u can be either one of the components (u x , u v ) or 
the velocity modulus. 

It turns out that inverse temporal statistics does not 
posses good scaling laws. Therefore, we refrain from giv- 
ing any quantitative statement while we concentrate on 
some qualitative properties showed by PDF's of temporal 
inverse events measured at the two probes, Pi n and p ou t- 
In Fig. 9 we plot, for the probe p in , various PDFs 
P ( T Su/ (t(Su))) 7 at varying Su, all re-scaled with their 
mean value (t(Su)). First, we notice that PDFs col- 
lapse very well, indicating the absence of intermittency 
effects. Second, between the peak and the exponential 



tails at large r, each probability density function ex- 
hibits, on a wide range of scales, a power law behavior 
P(t~8 u ) ~ (tsu)" 1 with an estimated exponent 7 *=s —I. 
On the other hand, PDFs measured on the probe outside 
the vortex, p out , show a very different qualitative trend 
(inset of Fig. 9). In particular, there is not any clear 
power law behavior. This indicates that very large exit 
events become less and less probable outside the probe, 
something which must have to do with the absence of 
very smooth fluctuations in the vortex background. 
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FIG. 9. Exit-time probability density functions P(ts u ) 
measured on the statistics of the probe p;„. The curves, cal- 
culated for different exit barriers Su, are normalized by their 
first moment. In the inset, the same for the probe p ou t- In 
each figure, the straight line is the power law behavior r^ 1 . 

Altough qualitative, the inverse-statistics properties 
allow to distinguish among different temporal statistical 
behaviors associated to different fluid regions. 



VI. CONCLUSION 

To summarize, we have studied inverse statistics mo- 
ments for signals with a more than smooth spectrum, 
i.e., signals which are diffcrcntiable and with non-trivial 
stochastic sub-leading fluctuations. We have shown that 
statistical velocity properties of 2D turbulent flows are 
not simply described in terms of the spectrum slope. 
From the exit-distance analysis it is possible to high- 
light a whole spectrum of more-than-differentiable fluctu- 
ations. These, being connected with laminar events, arc 
the strongest statistical signature of the large-scale co- 
herence. Experiments with different methods of remov- 
ing/pumping energy at large scales should be performed, 
to investigate the importance of large-scale structures in 
the inverse statistics of flows with different spectra. We 
have quantified laminar fluctuations also by using Second 
Differences, i.e., direct velocity increments subtracted of 
their linear differentiable behavior. We have found also 
in this case that more-than-differentiable fluctuations are 
not simply described by one single exponent. 
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As a final remark, we stress that inverse statistics pro- 
vide a completely new statistical indicator with respect 
to the standard direct statistics observables. We have 
shown that such method is necessary in all those cases 
where non-trivial fluctuations are sub-leading with re- 
spect to the differentiable contributions. Obviously, the 
same kind of analysis here reported can be extended to 
other temporal signals, applying the method to a broad 
class of natural phenomena. As an example, we just men- 
tion possible applications in situations common to clima- 
tology or meteorology where estimating the probability of 
persistent velocity configurations, or of any other dynam- 
ical variable, is relevant. As a perspective, an important 
generalization is the investigation of multi-dimensional 
signals by studying the statistics of <i-dimensional vol- 
umes between equispaced iso-surfaces. 
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TABLE I. Drag parameters p,/3 p , spectrum slope a, and 
the real space sub-leading scaling exponent, h = (a — l)/2 for 
the various numerical experiments. The value of each a has 
been obtained by a best fit in the region |fc| ~ (20 — 100) (see 
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slightly larger values of a are obtained. These discrepancies 
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DNS label 
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Pp 
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h 


A 





0.01 


3.26(6) 


1.13 


B 





0.10 


3.38(8) 


1.19 


C 





0.30 


3.74(8) 


1.37 


D 


2 


14.0 


3.26(6) 


1.13 
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